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ABSTRACT 


By examining the interaction between supernova CSN) ejecta and the 

various environments in which the explosive event might occur, we conclude 

tha*: aly a small fraction of the many SNs produce observable supernova 

remnants (.SNRs), This fraction, which is found to depend weakly upon the 

lower mass limit of the SN progenitors, and more strongly on the specific 

characteristics of the associated interstellar medium, decreases from 

approximately 15 percent near the galactic center to 10 percent at Rgal 

M.0 kpc and drops nearly to zero for Rgal >15 kpc. Generally, whether a 

SNR is detectable is determined by the density of the ambient interstellar 

medium in which it is embedded, We find that SNRs are only detectable above 

* 

some critical density (n^0,l cm” }, The presence of large, low density 
cavities around stellar associationdue to the combined effects of stellar 
winds and supernova shells strongly suggests that a large portion of the 
detectable SNRs must have runaway stars as their progenitors. These results 
explain the differences between the substantially larger SN rates in the 
Galaxy derived both from pulsar statistics and from observations of SN 
events in external galaxies, when compared to the substantially smaller SN 
rates derived from galactic SNR statistics, These results also explain the 
very large number of SNRs observed towards the galactic center in comparison 
to few SNRs found in the anti-center direction. 

Subject headings; nebulae; supemoyae remnants - stars; stellar statistics 
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I. INTRODUCTION 


The evolution of a SNR expanding within the "typical" interstellar 
medium (i.e., number density nM. cm , temperature TMO -10^K) has been 
studied with a great deal of detail and sophistication (for a recent review 
see Chevalier, 1977). In particular, detailed integrations of the basic 
shock equations have been performed (Chevalier 1974; Mansfield and Salpeter, 
1974), and the complex phenomena which appear during the transition from 
the adiabatic to the isothermal phase hava been addressed (Chevalier, 1975; 
Chevalier and Theys, 1975; Woodward, 1976; McCray et al. , 1975, etc.). 

Recognizing the inhomogeneous nature of the interstellar medium (ISM) 
considerable work has been carried out to study the effects of inhomo- 
geneities in the structure and evolution of the SN shock waves. For example, 
McKee and Cowie (1975), Sgro (1975), and Woodward (1976) have investigated 
the interaction of supernova shock waves with interstellar clouds , and 
McKee and Ostriker (1977) have examined the effects of SN explosions on 
a cloudy interstellar medium. 

In all the above work, the basic shock is always assumed to propagate 
within a typical (as defined earlier) interstellar medium, and until quite 
recently, no attempt had been made to study the evolution of a SN shock 
wave expanding into media representing the various possible environments of 
supernova progenitors within the Galaxy. For example, it is well known that 
most (if not all) stars are born in groups (clusters or associations). 

The precursor of the stellar association is a dense molecular cloud. The 
first supernova from the stars in the group is set off near, perhaps 
inside dense, cold gaseous media. The evolution and long term detectability 
of such an event (recently studied by Wheeler et al. , 1980; and Shull, 1980) 
is obviously very different from that of the canonical supernova remnant 
that propagating in the typical interstellar medium . 
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Supernova* propagating is deua* molecular clouds will be very rare, 
however, since the first on* nay disrupt on* of the nailer clouds 
(Wheeler et al . 1980) , whereas for th* more massive clouds, the combined 
affect of stellar winds and the earliest supernova* (Bruhweiler eit al . , 

1980) creates an expanding, hot, low density cavity, within which subsequent 
supernova shell* will expand for times up to tens of thousands of years. 

It is this scenario, in fact, which is by far the most common encounters' 1 
by supernova sheila. As we shall discuss in ^ V, however, since most 
massive stars occur in binary systems, when the primary member of the 
system becomes supernova, it may impart a large velocity to the secondary 
star by means of the slingshot effect (Blaauw, 1964) . Many of these stars 
will in fact overtake the supershell and remain within the confin.s cf 
the galactic disk until they become supernovae. These stars are the ones 
that produce the typical supernova remnants so extensively discussed in 
the scientific literature. 

In this paper, we will study the evolution of supernova shells into 
these very dissimilar media. Because of the dramatic effects of the pre- 
supemova environment on the evolution of the remnant, we have confined 
our theoretical discussion on the. simplest possible description. Thus, 
following Spiczer 0-978) , we characterize the supernova shock by three 
phases, namely: 

(a) initial free expansion of the supernova material, 

Cb) Intermediate adiabatic (or Sedov) expansion, and 
Cc) late isothermal expansion. 
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Although the structure of the pre-supernova environment can, in general, be 

fairly complex, it is usually made up by a combination of dense clouds, 

diffuse, warm medium, and rarefied, hot cavities. To that extent, we shall 

study the evolution of supernova remnants into these three types of media, 

which should span conditions for nearly all real cases. Thus, we will 

consider the remnant evolution into 

4-3 <2 

(1) a dense (n ^10 cm ), cold (T ~ 10 K) molecular cloud 

(2) the typical interstellar medium (n ''d cm’ 3 , T MO^-IO^K) , and 

-2 -3 5 

(3) a bubble or supernova cavity (n ''dG pm , T ^5 x 10 K). 

II. TYPES OF SUPERNOVA REMNANTS 


Before describing each particular type of SNR let us introduce the 
relations which allow us to compute the several SN phases for each particular 
environment. As stated earlier, we shall follow the formalism of Spitzer 
(1978), and Gorenstein and Tucker (1976). 

In the early phase, (a), 

* shock wave will travel just ahead of the ejected shell, with a velo- 
'city Vg. The shock will heat matter to a temperature 

T s = 3^ M Vg/16k ^I.3clO" 9 Vg (l) 


(Landau and Lifshitz, 1959), where is the mass of the hydrogen atom, y 
is the mean molecular weight and k is the Boltzmann constant. This phase 
ends when the swept-up interstellar material equals the mass of the ejecta, 
i.e. , when 



where M^ is the mass of the ejecta and p the density of the ambient inter- 
stellar medium. The elapsed time, t, between the supernova event and the end 
of phase (a) is given by 
R S 

t * t r o) 
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Phase (b) , the intermediate non-radiative expansion phase, can be computed 
by means of the Sedov solution. The temperature immediately behind the 
shock is given by 
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where E is the kinetic energy of the ejecta j 
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and t can again be obtained from (3). 

Phase {b) ends when T falls below M0 6 K, for n 'v 1 cnf^ since radiative cooling 
then becomes important, which brings the onset of phase (c), the late Isothermal 
expansion. This phase can be represented by the snowplow model, where 
conservation of momentum applies. Here the shell velocity is given by 



where and are, respectively, the shell mass and velocity at the end 

of phase (b). At this phase, most of the swept-up material is in a cool, dense 

shell although there may be sore thermal X-radiation from the hot low density gas 

interior to the shell. This low density interior gas has a long cooling time. 

We now discuss the three example SNRs and how differently they appear 
* 

in the several expansion phases. In all of our calculations we shall assume 

that during the SN event, 4 M of the stellar material are ejected with a 

© 

velocity of 5000 km sec ^ , and thus, the total kinetic energy of the ejecta 
51 

is 10 erg. These parameters are a realistic representation of a type II 
SN which would be expected to occur in an 0B association. 


A 


A. A SH in a Dense Molecular Cloud 


While molecular clouds found in nature have a large range of size, mass 

and density, a typical molecular cloud can be approximated as a sphere with 

4 -3 

uniform density n^ ^10 cm having a diameter D c ^<5 pc (i.e., or cloud mass 
M c ^3. 5 x 10 4 ), and with a temperature T < 10 2 (Burton, 1976). 

Wheeler et al., (1980) And Shull Q.980) have modeled the effects of 
a SN exploding within a molecular cloud. The characteristics of the SNR 
at the end of each phase using their results are summarized in Table 1. . 

We will subsequently call these SNRs the molecular cloud SNRs . 

Because of the very high density, the phases occur rapidly for the 
molecular cloud SNRs. Wheeler et al . (1980) indeed suggest that the adiabatic 
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phase may not exist if n„ > 10 cm . Whether this happens depends on n as well as on 

2 

M e j. The molecular cloud SNR represents the conditions within a young mole- 
cular cloud when the first, most massive stars become SNs. These conditions 
wou!a be very short-lived as the SNRs would push the molecular cloud away 
from the remaining massive stars within the association and cause rapid 
cloud fragmentation (Elmegreen, 1979). Consequently, this picture will apply 
to a relatively small minority of the SNs which occur in the Galaxy. More- 
over, at no point of its evolution does the resulting shell remotely resemble the 
familiar observed supernova remnant (SNR). In particular, during the earlier S'.ages 

of evolution, optical and x-ray observations are useless as a means of 
detection because the cloud is optically thick to those wavelengths. If the 
SN were to produce a gamma ray pulsar as the stellar remnant, it would be 
observable. However, since only a minority of the pulsars are known to emit 
gamma rays, this is not an effective means to search for SNs within dense 
clouds. The SNR should be detectable by means of the infrared emission from 
the heated grains inside and outside of the cloud (Wheeler ejt aK , 1980; Shull 1980; 
Silk and Burke, 1974). Even this technique is not fool-proof, though, since it may be 
very difficult to differentiate the cloud-embedded source of infrared emission 
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as a SNR rath*? than a recently formed OB association, 
i. A SN Within a Hot, Rarefied Cavity 

Most SNs which occur in the older OB associations will first expand into 
the hot highly-evacuated volume produced by the combined effects of stellar winds 
and earlier SNe (Bruhweiler et al. , 1980). The ejecta expand freely until 
enough gas is encountered to form a shock. The ejecta do not encounter • 
significant gas until very late in the evolution of the SNR. Indeed, for a 
few thousands of years the SNR is hot, dilute gas expanding without bound. 

Such a gas is very difficult to observe (the question of detectability will 
be addressed in<f III. 

As an illustration, we compute the evolution of a SNR contained within a 
superbubble with radius R * 105 pc. From the model calculations of Bruhweiler 
et al . , such a shell would exist around a typical OB association after a few 
million years. The characteristics of the SNR within a superbubble are sum- 
marized in Table I at the end of each phase. We shall call these SNRs the 
hot cavity SNRs . 

C. A SN Surrounded by the Undisturbed ISM 

In some instances, an intermediate mass star may become a SN outside 
molecular clouds and also outside the hot cavities surrounding OB associations. 

The resulting SNR which interacts with the previously undisturbed interstellar 
medium of the galactic disk is the canonical SNR that has been described in f I. 
The characteristics of the canonical SNR are summarized in Table 1 . 

The canonical SNR has several intriguing differences in properties when 
compared to the other two types of SNRs. The total evolutionary lifetime is 
substantially longer than the lifetimes of either the molecular cloud SNR or 

the hot cavity SNR. The kinetics of a molecular clouu SNR are quickly transferred 

4 

to the very massive moleculur cloud and the expanding shed 1 stalls within 9 x 10 
years for n *'» 10^ The hoc cavity iNR, on the other h-.r>-‘ , expands rapidly 
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until it encounters the outer, slowly moving supershell of neutral gas 
(as observed by Hailes, 1979). By then the SH shock is highly diluted, 
and consequently the ejecta are quickly decelerated. The slow accretion 
of material by the classical SNR in the undisturbed ISM extends the 
lifetime by more than twenty- fold compared to either alternate example. 




The EM observability criterion (EM £ 50 cm' 6 pc) is a function of time. 
For the Sedov phase, we find 

■o l »•« E 5l' 1/5 h 2 ' 9 c" 3 (8) 

where EM > 50 cm' 6 pc. Taking E^. -l.yO (see part C of Table 1) we find 

that n £ 1.4 era J ie required. even lor the SNR with 
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ui ^ s cm pc, we need > G.J cm 

For an SNK in the isotherms 1 phase , f.v- -oM^v.sb} ’ vcy estimates 


are different. Large comp 


' ' ey estimate! 
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111. OBSERVABLE SWRs 

We discuss here the type of interstellar environment required to 
produce an observable SNR. By "observable” we w*->tn an SNR which may be 
seen in a) visible light, b) radio waves, or c)X-rays. 
a) visible light SNRs 

The majority of visible light SNRs are believed to be in the Sedov (adiabatic) 

phase. Exceptions may be very "old" SNRs like the Monoceros Loop, although McKee 

eui Cowl* (1975) have suggested that even those sTe in the Sedov phase, or ver, 

"young" SNRs like A and the Crab Nebula. A SNR can be easily detected at visual 
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wavelengths if the emission measure, EM * n g L, exceeds 50 cm pc where n # Is the 

electron density to the SN shell and L is the shell thickness. Careful observations 
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will aid in detecting an SSTi with EM 20 -50 cm pc for T v 10 K but very special 

tedious techniques are needed to detect a SNR with ©I o# 5 cm 6 pc (for example e 
large Fabry-Perot etalcn, etc.). 

The emission measure can be expressed as (assuming, roughly, an average path 
length of R fi /12 in the Sedov phase) 

EM = t K Til (7) 

3 S o 

where n o is the ambient ISM density and Rg is the radius of the shock 
front. Equation (7) is actually an over-estimate because the density 
drops off rapidly behind the shock front (Spitzer 1978; chevalier, 1974). 






find that compression ratios as high a* 'v 70 nay b« rssohsd if tha ambient 
magnetic field is weak or as low as * 20 whoa tha aagnstio fiald oaanot bt 
ignored (1 >3 X 10 a). Tha appropriate pathlangths is tha characteristic 
oooling length (McCray at § 1 , 1975). For a 100 ka s’* 1 shook with no 
anbient magnetic fiald (tha most favorable ease) tha aaission measure would 
exceed 50 cm’ 6 pc if » o £ 0.9 cn“ 3 . Correspondingly , for EM £ 5 oaf® pc 
B o as" 3 . Since the aabiant densities are lower, it is wore favorable 

to o b se r v e an optical SNR during its la otho r aa l phase than during its 


adiabatic phase. For an "easy" detection, it 


that the shock Bust 


collide with a dense cloud (locally, the intercloud medium has a density 
of * 0.15 cm* 3 , Falgarone and Laqueux, 1979). 

The filling factor of clouds is quite mall (f 1-10% see McCray 
and Snow 1979). Moreove r, the number of clouds drops off rapidly with 
height above the galactic plane (the cloud scale height is probably <v 1/2 
of the male height of the diffuse intercloud medium; see Falgarone and 
Laqueux, 1979). Hence, few optically observable SNR's are expected at 
large heights above the galactic plane. The Cygnus Loop say be observable 
because it has collided with a neutral cloud or clouds. Thsse clouds 
cannot be the sane type of neutral clouds that product H I absorption 
profiles, as tha H Z absorption clouds ire found (Radhakrishnen et al., 

1972) within 300 pc of the galactic plane, 
b) Radio observable SNRs 

Radio SMRs are such a ore numerous than either optical or X-ray reanants . 
About 130 radio SNRs have been observed in the galaxy (Clark and Caswell 
1976), whereas only 30 of tha radio SNRs are detected by any optical emission 
(van den Bergh, 1975). It is easier to dettet a SNR at radio wavelengths 
both because lower ambient densities are necessary for radio detection and 
because dust does not absorb radio photons. Zt is usually assumed that 
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r - 1 ' 1 
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van 4ar Laan's (1962) theory applies to the older SNRs. However, the 
statistical investigation by Clark and Caswell (1976) confirmed that for 
the majority of SHRs the Sedov solution properly describes the value of the 
diameter D with time. Clark and Caswell derived an average <f7n Q > * 5 X 10 
ergs cm 3 where E is the initial supernova energy end n Q Is again the ambient 
density. If E ■ 10 erg, the averaged n Q implied by their analysis is 0.2 an . 
More recently, Caswell and Lerehe (1979) refined the t - D relation to include a t 
dependence. (E la the radio surface brightness.) The derived scale height of rad 
SHRs la 200 pc. The lulled frequency of SN producing radio SNRs within the Galax 
if f 1/80 yr, * The lower values of the ambient density implied by the radio 
observations confirm that SNRs are mote easily detected In the radio. 

There is also information on the galactic distribution of SNRs. Tlovaisky 
and Lequeux (1972) find that the distribution of radio SNRs closely follows 
the radial distribution of the non-thermal background radio emission. At 
t kpc from the galactic center the radio SNRs are three times as abundant 

at 10 kpc. Beyond 15 kpc there are very few radio SNRs. 
c) X-ray observable SNRs 

Very few X-ray SNRs have been observed (for a recent revi ew of’ the 
X-ray SNRs sec Clark and Culhane, 1976) perhaps being due in part to the 
limited sensitivity of complete X-ray surveys. The E/n value derived for 
the X-ray SNRs support the average value obtained from radio observations 
(Gomnstsin, Harnden and Tucker, 1974; Clark and Caswell 1*176) although the 
two studies are vastly different. The thermal X-ray flux in the .1 1 00 keV 
region of the spectrum for SNRs is given bv Goreustein, Harnden am Tucker 

Using the above expret f r > v m,d the 
X-rajfc measured for Pup A, the Cygnus Loop and the Vela ' '■ . >«f <;,-n«-lud* than an 


as L v "t 5 X 10“ 3 R 3 erg 
x o pc 


X-ray SNR is not observable if the ambient density is *: 

-3 


than 'v 0.1 cm 


On the basis of r« i? HRAO - 7 ("v < 
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for X-ray emission of radio and optical SNRs (Knox Long, private communication^ 

initial ambient densities below 0.05 cm“^ would not produce detectable 

X-ray fluxes. A value of the critical density of 0.1 cm i^vlred to observe 

an X-ray remnant is consistent with the "Einstein” results. From the above 

discussion we conclude that if the ambient interstellar medium density n Q 

exceeds a critical value n c , the SNR is observable. It would most probably 

be observed as a radio SNR, but if n^ is appreciably larger and the SNR is not 

too distant, it may be an optical or an X-ray observable SNR. We adopt the 

value n *0.1 cm~^ with the awareness that this is a realistic estimate really 
c 

for the radio SNRs. In any case, we find that the study which follows would 
not change appreciable If n were to change by a factor of 2 to 3. For optical 

C 
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SNRs in the Sedov phase, we find 'i* 1.0 cm but this value Is probably 

unimportant since it is much easier to detect the radio SNR. 

To simplify our analysis, we assume that al.l SNRs located within a 

medium with ambient density exceeding n will be observable, but all SNRs 

c 

located where n Q <n fi are not observable. For example, SNRs with diameters 
approaching the cloud scale height would be expected to be brighter on 
their edge nearest to the plane (Clark and Stephenson, 1977; Caswell and 
Lerche, 1979) and therefore, would be observable even though a portion of 
the SNR is located in a medium with density less than n . The largest 
observable SNRs have diameters < 50 pc which is less than the cloud scale 
height, so our results will not be significantly changed by this effect. 
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IV. STATISTICS or SUPERNOVAE 


The average expansion rate for OR associations Is ~ 5 km s’*. Since 
the radius of a typical supercavity at the Sun's distance from the galactic 
center is 250 pc, it Cakes about 5 x 10 7 years for an association member to 
overtake the supershell. During this time all stars more massive than 7 
will have completed their evolution, so that if a lower mass limit of 8 M. 

V 

• ’¥* "'jail & ' •' *$il •: 1 ^ £'< ’ i - 1 '» ,t > V* . 

for stars that produce SNs Is adopted, we find that very few stars are 
capable Of Escaping the cavity before they become supernovae (the numbers 
are higher for Closer distances to the center— see below — but still small). 


If the mass limit of stars producing SNs is extended down to 4 M . then 


the majority of SN explosions that produce SNRs in the solar neighborhood 


would originate in 4-8 M . A significant fraction of the progenitors 


would escape the low density cavity before exploding. As there is an 


uncertainty of the mass limit of stars which produce SNs we will consider 
two limits: 8 M and 4 M $4»t£ln part (A) below, we use 8 M while in part 


(B), where we also include the effects of the gravitational field from the 


galactic disk, we adopt 4 M* as the lower mass limit. By these two examples 


we represent the upper limits and lower limits to supernova occurrences 


within OB associations 


Early type stars have a high incidence of blnarie 


early B stars (Abt and Levy, 1978) shows that about hul 


multiple systems 


Rlaauw (1954) suggests chat when a SN occurs in such a binary system 

the companion can become a runaway star. In Ms classical work, 19 

runaways are identified out of which the latest spectral type is B3 with 

an assigned mass of ^ 1C , A later study of 304 0 stars (Cruz et al_,, 1974) 

- concludes that at least j. r ' percent of all 0 scats ate ;ur,;nv stars. Using 

» * 
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the above evidence, and assuming that all blnarlea produce runaways, it 
follow* that one-third of ell SNa are froa runaway progenitor a . It la 
the runaway star* that have an opportunity to eacape the auperbubblea 
and In turn produce observable SNR*. 

However, not all runaway progenitors produce observable SNR*. In 
order to estimate the fraction that produce observable SNRs, we need to 
estimate the following: 

I) the Initial mass function (IMF), 

II) the total evolutionary lifetime, T, for stars with different masses, 

III) the average peculiar velocity for runaway stars, Vp, which when used 
with the total evolutionary times will allow us to estimate the 
total distance, d, that a runaway star would travel froa the OB 
association, and 

iv) the effective critical scale height, H^, for observable SNRs. 

Ve now discuss each of these variables that will Influence the estimate 
of observable SNrs. 

The IMF for the massive stars is uncertain at best, expeclally for 
the 0 stars. Ostriker et al ., (1974) used the observational data of Richstone and 
Davidson (1972) to derive an IMF for 0 stars. Their mass function 

predicts significantly more massive 0 stars than that predicted by an 
extrapolation of the IMF deduced for the mid to late-B stars by McCuskey (1966). 
Theoretical evolutionary calculations, when compared with eclipsing binary 
data (Stothers, 1972), indicate a different mass versus spectral type 
relationship than used by Ostriker et al. Bruhweller (1980) has re-analyzed the Rich- 
stone and Davidson (1972) data and has determined the masses for MK 
standards based upon the work of Stothers (1972). Both the data of Richstone and David 
son and of McCuskey (1966) can be represented by an IMF which relates 
the total number of stars as: 
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an OB association before it becomes a supernova: 


where T. and T* are the total evolutionary times of the primary and the 


secondary in a binary system and At,. is the difference. The average 


mass of the primary, Mj, can be expressed in terms of the mass of the 
secondary, M^, by the relation: 

Ai /» \“2»21 i. L ft 


The higher mass limit for the IMP is not critical. In equation (11) 


ve assume that the primary has a random mass distribution described by 


the IMF in equation (9). In Table 2, we present the resultant mass 
(which is a mean value) of the primary, Mj, computed from equation (11) 
and the timescale, At,. » T. - T,, which enters the expression in (10)*. 


FOOTNOTE 


*In our discussion we have ignored the effects of mass loss and mass exchange 
in the evolution of 0 and B stars. These processes affect A x 12 , which in 
turn affects the fraction of runaway stars that escape from the supercavity. 

Mass loss is expected to lengthen (by about 102) the evolutionary lifetimes 
of the more massive stars (Chios!, Nasi and Srecnivasan, 1978). However, the 
domain of extensive mass loss is limited to bU or earlier stars (Snow and 
Morton, 1976), i.e., stars with initial masses > 17 M ft . From our calculations 
these stars become supemovae either inside the supercavity or outside the 
galactic plane. Thus, the proposed increase in evolutionary timescales has a 



negligible effect. 

A potentially acre significant effect is due to the mass distribution of 
the binaries. For the sake of simplicity, we have assumed a random distribution 

of mass ratios among binaries. This assumption gives Mj/M^ > 0.73 for masses 
of interest in this work. On the other hand, if close binaries are formed 

as bifurcation products, then the mass ratio should be of order unity. 
Observatlonally determined values for this ratio range from 0.35 (Stone, 1979) 
to about 1 (Helntse, 1973). In view of these uncertainties, jur assumption 
of Mj/M^ ^ 0.73 is reasonable. If the true ratio were smaller, A 
would be larger than our tabulated values, whereas the bifurcation hypothesis 
leads to smaller than our tabulated values. In 

any event, the effects are not overwhelming, and our procedure is justifiable 
at present. A final uncertainty is introduced by our having ignored the 
effect of mass transfer incomputing A Due to the speeding up of the 

evolution beyond core H~exhaustion stage, the effects of mass dumping by the 
primary is expected to be minimal when compared to the shortening of the 
evolutionary timescale of the secondary. This would lead to generally smaller 

^ T 12 * Since in otder to produce a runaway the mass ejected must exceed 

the mass of the companion, mass transfer could modify the statistics of close 
binaries . 


The effective critical scale height, H c , for observable SNKs is 
determined by the scale height of the gas. In Figure 1 we shou the 
structure of the supercavity produced at three different distances from 
the galactic center, R ^ ■ 5, 10, and 20 kpc. These are structures based 
upon the model we calculated (Bruhweiler et al.» 1980). Ue assumed an 
exponential density distribution 


n ■ n e 
o 


-*/H 


( 12 ) 


with H • 70, 150, 500 pc and n Q ■ 3, 1, 0.1 era ^ for R^^ ■ 5, 10, 20 kpc, 
respectively. The assumed densities are appropriate for the H I medium 
as determined in Paul, Casse' and Cesarsky (1976) and the scale height is 
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from Kerr (1969). Exponential distributions were found by Celnik, Fohlfa 

and Braunefurth (1979) for large distances away from the galactic plane. 

They give detailed formulae for the scale heights, H, although the values 

of are harder to determine from their work. Assuming that 5 percent 

of the total mass of the galaxy Is In gaseous form, we find - using the 
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work of Celnik, Rohlfs and Braunsfurth - that n Q ^ 6, 0.7 and 4 x 10 cm 
for ■ 5, 10 And 20 kpc respectively. Bohlin et el. (1978) on the other 
hand, find that n Q 0.9 for the H I medium and ^ 1.2 for the H I + medium 

(in the solar neighborhood). 

Even though there are large uncertainties In the parameters of the ISM 

gas, the structure of the bubble-SN cavity is not affected very much. This 

is so because the radius of the bubble is only weakly proportional to the 

ambient density (“ n ”*^) and the radius of the SN produced shell is a n 

o o 

near the plane. For the case where the uncertainties are large (R^ “ 20 kpc) 

-3 

we find that even the higher density U q m 0.1 cm produces such large super- 

cavities that essentially no runaway stars escape. Therefore, this result 

would still be true if we choose n to be smaller. 

o 


A. Supernova Statistics Derived with 8 Lower Limit for SN Progenitors, 
and Ignoring Gravitational Effects of the Disk. 


We now compute the percentage of all SN progenitors that produce observable 
SNRs. In this case we assume 8 to be the lower limit for SN progenitors, 
and ignore the gravitational effects of the dick. The gravitational force 
of the disk tends to restore gas to the galactic plane and in the ? -direction 
decreases the size of the supercavities. Hence, the percentages found in 
this case will be a lower limit to the actual percentage. 


1G 


«.< m* *» - ■ •*» 


In Figure 2 we present the fraction of runaway etare in mate intervale 

AM • 1 Mq that escape the eupercavity to produce observable SNRs. The 

results change for different ratios of H^V^; hence, we show the curves for 

a range of H /V . For example, if we assume V ■ 50 km and we want to 
c p p 

appropriate curve for ■ 5 kpc, from Figure 1, H £ ■ 240 pc and therefore 
H /V ■4.8, Similarly, for R-- ■ 10 kpc H ■ 350 pc and H /V ■ 7. 

C vAX C C p 

We now illustrate the application of this figure by presenting numerical 

estimates for the fraction of stars in a 1 mass interval that escape the 

“z/H 

supercavity and produce an observable remnant. We chose n - n Q e as in 
paper 1 at the two galactic distances R Q|k ^ ■ 5 and 10 kpc. The critical 

_3 

density, n £ ■ O.i cm , makes the total fraction of SN producing observable 
remnants at R^ ■ 20 kpc equal to zero. Decreasing the critical density to 

_3 

0.05 cm would predict a total fraction at R^^ * 20 kpc of less than 1 percent. 

Table 3 presents the lower limit on the percentage of SNs that produce 
observable SNRs. In Column 1, rows 2-10, we Increment the mass range in bins 
of 1 Mq. However, row 1 has a bln of 17-70 Mq . In column 2, we present 
the fraction of all SN progenitors that are within the mass bln. In columns 
3 and 4, the fraction of all SN progenitors, which are within the mass bin, 
is given for those that produce observable SNRs at 5 and 10 kpc. Column 5, 
which would list the fraction of SN progenitors that produce SNRs at ■ 20 kpc, 

is empty to emphasize that no SNRs would be produced in the ambient interstellar 
gas at 20 kpc . In row 11 we add the incremental percentages to find the 
total percentage of runaway SN progenitors that produce observable SNRs. To 
this total, we must add n row 12 the few SN progenitors in the low mass range 
that survive long enough to escape the supercavity even at the association 
expansion velocity of 5 km s By comparison, we have assumed that the 

runaway progenitors have ■ 50 km s~*. We see that no slow moving stars 
escape the supercavity at ■ 10 kpc while only 4,3 percent of the SN 
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progenitors are slow moving stars which escape the supercavity of ■ 5 kpe. 

Even if a SN explodes outside the supercavity formed by the parent 
association, such a SN might not be in the ambient Interstellar gas. Rather 
the SN may find itself in another supercavity and hence, it would not form 
a visible SNR. The fraction of the ISM occupied by these supercavltles is 
hard to estimate. The 0 VI gas (McCray and Snow, 1979) has a filling factor 
estimated to be twenty percent whereas the hot gas responsible for the X-ray 
background (Krmishaar 1977) has a filling factor estimated to be 50 percent. 

We find that about 30 percent of the ISM is occupied by these superbubbles 
produced by OB associations. Bowever, this is likely to be a lower limit 
since we do not include In our model the (older) B associations. With 30 
percent of the ISM assumed to be in the hot phase, wc find the percentages 
listed in line B. Note that we have decreased the percentage of the runaway 
progenitors by 1/3 and not the percentage of the slow-moving progenitors, 
as the latter will be just beyond the superbubble and would be very unlikely 
within another supercavity. 

B . Supernova Statistics Derived with 4 Lower Limit for the SN Progenitors, 

and Including Gravitational Effects of the Disk. 

The runaway progenitors and the gas are subjected to a gravitational 
restoring force towards the disk which we have ignored thus far. Close to 
the galactic plane the gravitational force lav can be approximated by 
z ■ -kz with the resultant motion being that of an undamped harmonic 
oscillator . 

We now calculate the critical angle, 9^, at which runaway progenitors 
would produce SNs at a height above the plane for a selected mass range 
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as represented by At 2 j ■ T 2~ T r In tera * the z-dlrection grevltetlonel 
force, we cen express z end H c s 

I - V p sin 0 H - kfk (13) 

end 

H c - A sin (fk At 21 ). (14) 

We cen solve for 0 H as: 

sin6 H -Si5 |,in (fk At 21 )| -1 (15) 

Within the ness range represented by A T^ . all runaway progenitors thst 
are ejected from the galactic plane at ejection angles less than 0^ will 
produce SNs at heights less than H c above the galactic plane. The fraction 
of runaway progenitors that produce observable SNRs for a given mass range 
AM 2 is then: 


H 

(16) 

The fraction f CM _ can be evaluated by substitution of equation (16) into 
equation (15). 

We are aware that large uncertainties are inherent in the gravitational 
force law for the Galaxy. Although errors in the force could be quite large 
at high latitudinal distances, we are mostly concerned with z-distances less 
than 300 pc. The z - -kz approximation has estimated errors of less than 

20 percent (see paper 1) in the solar neighborhood (R Cal * 10 kpc). However, 
at 5 and 20 kpc the force law is much more uncertain. 

Based upon the curves published by Schmidt (1956), we adopted in 
Bruhveller at al. (1980) linearized force laws out to 300 pc in z with the 
values of k being 6.06 x 10 - **, 2.58 x 10 and 3.65 x 10 ** sec ^ for 


\ H 

2 ’J -e a <v , fit 2i 


(AM 2 )) cos 0 d 6 


SNR 


sin 0 


At 21 (AM 2 )) 4 


a* 
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*G*1 " 5 * 10 and 20 kpc respectively. 

Two other uncertainties are tha values for B. and V.. The derived 

* « P 

values for H depend upon the assumed values of n^ end the assumed density 

''*.7 ' 

lew (Equation 12). The average velocity, V p , for runaway progenitors is 
also very uncertain. Several runaway stars are known to have >_ 100 las s 
&*»uw, 1964s Stone 1979). The average runaway velocity we use is based on the 
velocities of 19 runaways presented by Blaauw. Various selection effects 
ere known that strongly suggest the lower velocity runaways may not be 
recognisable as such. While a reasonable velue for may be on the order 
of 50 km s“ l (or lower), we feel constrained to express H c and V p in 
terms of one variable, namely, their ratio, H c /V p . 

We now show in Figures 3 and 4 curves similar to those in Figure 2, 
but with the gravitational restoring force included. Within the mass intervals, 
A M - 1 Mg, , we plot the fraction of runaway progenitors that escape the 
supercavlty and produce observable SNKs, Figure 3 is for w 5 kpc and 
Figure 4 is for - 10 kpc . We present the results for the ratios of 

H /V - 1.0 and 2.0. 

C p 

The percentages of SNs that would produce observable SNRs are summarized 
in Table 4. As in part (A) we assume that 1/3 of all SNs are from runaway 
stars and that 1/3 of the interstellar medium is occupied by supercavities. 

The percentages listed in Table 4 are for an evolved supercavity with the 
internal SN from low-mass progenitors being at large distances from the shell 
(labeled in Figure 1 as SN2). The SNR from the more massive stars have 
either dissipated or have been overtaken by the expanding supercavity around 
the old association. If stars appreciably less than 8 Mq produce SNs, then 
it is possible that these older supercavities would be ringed with SNRs from 
the slow moving progenitors. For example, if the expansion velocity were 
5 km s' 1 for the association, then 5 Mq stars would diffuse 500 pc, which 


?0 



| 


In the solar neighborhood li twice the shell redine, clone the plane, by the 
tine they becone She. 

Ve have also computed c la liar models for young eupercavltles where only 
stars with spectral type BO or earlier have become SN (the radius is shown 
in Figure 1 as "SHI"). These younger eupercavltles would be surrounded by 
SNRs from the relatively more massive progenitors (M > 11 Mg ) . A possible 
example of such e system may be the dm Nebula. It is a roughly spherical 
cavity with a 125 pc rad in a. Two known SNRs in close proximity are Vale X-l 
end Fuppis A. 

Tim computed estimates of SNs producing SNRs for the relative younger 
and older supercavities are very similar. The total percentages, summarized 
in Table 4, are 29.8 percent and 23.1 percent for R^ ■ 5 kpc and 10 kpc 
respectively. About one half of the SNRs are from the alow moving (5 km s”*) 
stars. 


Ve also computed the percentages of SNs creating observable SNRs at 
W 20 kpc and find that less than one percent of the SNs would yield 
observable SNRs. 
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V. DISCUSSION AMD CONCLUSIONS 

Several important results are derivable froa Tables ^ 

The presence of supercavities drastically changes the mas distribution 
of SM progenitors than produce tha classically detectable SNRa. 

A greater percentage of higher-mass stars escape the supercavity 
at snail R^ primarily because the eises of the supercavities are smaller 
toward the galactic center. The nodal prediction that acre SN produce 
detectable SNKs at wall R c#1 is in qualitative agreement with the observed 
distribution (van dor Bergh, 1978; Clark and Caswell, 1976; Ilovaisky and 

_3 

Lequeux, 1972). Had we chosen n Q to be 6 cn — which stay be more 

appropriate for the dense inner eras— we would have predicted 2.2 tines more SNRs 
at 5 kpc than at 10 kpc, Due to tha uncertainties involved we consider 
this to be satisfactory agreement .with the observations. This also implies 
that our assumption that the rate of SN outbursts throughout the galactic 
plane us nearly uniform way be close to reality. The distribution of 
observable SNRa is only determined by the presence end else of the super- 
cavities end by the fraction of "runaway” SN progenitors that escape before becoming 
SNs. 

The exact fraction of SNs that produce observable SNRa cannot bt 

estimated accurately with the present data. Upon a variety of assumptions 

and parameters hopefully encompassing the real situations, we find that 

this fraction ranges from approximately 10 percent to 35 percent depending 

on the with 5 percent variations In either direction being reasonable. 

The Interstellar gaa densities effect the above percentages somewhat. 

However, the fraction of SNs producing observable SNRs more than doubles 

when the lower mass cutoff decreases from 8 HL to 4 m_. Moreover, if 4 

O O Q 

stars ara progenitors, the slowly moving stars which diffuse at the association 
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velocity of expansion (5 ka s ) would make an equal contribution to the 

total SNR percentage. 

Since our model indicates that only one to three out of every ten SN 
produce an observable SNR, the dlsagreesMnt between rates deduced from observable 
SNRs on the one hand and SN end pulsars on the other hand can be understood. 

Tamaann (1974) suggests e mean interval between SN,t sn *v 30 years, this leads 
to a SNR production rate of one every 100 to 300 yr. On the other hand pulsar 
statistics (Taylor and Manchester, 1978) imply t $n % 10 years. This leads to a 
SNR production rate of one every 30 to 100 years. Since the mean T gNR derived 

by Caswell and Lerche (1979) is ^ 80 years, there is a weak support for the 
higher SN rate. However, regardless of whether one SN occurs every 10 or 30 
years, we can now understand how both these numbers are lower than the 80 year 
Interval for the production of detectable SNR derived from radio data. 

Recently, Higdon and Llngenf alter (1980), have proposed an alternative to 
our point of view. They propose that if a hot (^10 6 K), tenuous gas fills 90X 
of interstellar space, the observed number and surface brightness distribution of 
galactic remnants implies e SN rate oi once every n. 30 yr. However, even though 
in this way the statistics may equally be reconciled, we feel that their filling 
factor is excessively large, although not inconsistent with the 0 VI date (Jenkins 
and Meloy 1974). 0ur Bech * ni8a would operate in 

eny case, further reducing the number cf observable SNRs. Bec.u.e of this, the 
90X filling factor is not Justified. The results of Tables 3 end 4 would 
indicate that filling factors as large as 70X could be tolerated but not 
appreciably higher than that. 
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We auet emphasise that our discussion has been confined to type 11 
supernovas. Type I supernovas, which probably originate Iron low-mass 
popular on XX stars, share alrwst none of the considerations addressed in 
this paper, and they should be investigated separately. However, since mbs 
statistical studies of SN do not dif CerentUtS the SN types, they cannot 
be compared with our results wltliout previously estimating the type X SN 
contribution. If one, however, examines the known SN of the last mlllenlum, 
all of which left observable remnants, the question of differentiating between 
SN of type X end SN of type XI becoaea important. 

Wt conclude that the aupercavlty concept— that of a hot, low-density 
gaa around stellar associations created by stellar winds and SNs — provides 
a most effective scenario for understanding the general structure of the 
interstellar medium. Moreover, we bviieve that supercavities have already 
been detected observatlonally as the HI supershells (Heiies 1979) and as 
the X-ray superbubble in Cygflus (Ca sh et al 1980). This provider a 
compelling argument that many SN do not occur in environments conducive to 
detectable SNRa . 


24 


,,1,1 f iwih'iii mi iw ^iWi"i | !w:<ii! i wsii wwik»i'tii‘ i,iir.-pn;Mrji ;::!:ii;iiuiilf ryiurnnir’y rfl)ll^T*til^|i 'I ' i ^ r swr, mwim-wniiwinni «•■.». w '•""■h 


1000 


H-150 pc 
n c =1 cm-3 


10 kpc 


n = n 0 e' 


1000 1 


H= 70 pc 
n 0 = 3 cm* 3 


R«„i =5 k P c 


•j 

a> 500 

N 

400 


0 100 200 300 400 500 600 700 800 9001000 
R (pc) 


a 500 


SN2 

£SN1 


-0.05 




0 100 200 300 400 500 600 700 800 

R (pc) 


Figure 1. Geometry of the supercavity stnictiues at 8 *m 

stages of evolution are depicted. The cum ^^ed B is the tout Jh ^hTewlutionary time when aO 28 massive O stars have 

S5S SNs 5 ** HaswreryiSo'mM^B^rt^^w^maa^ ^Je^ter than 8 M©. These too become SN, and drive the supercavity to Ore 

SN2 size. 





LOG OF PERCENTAGE FOR MASS INTERVALS OF 1M 







LOG OF PERCENTAGE FOR MASS INTERVALS OF 1M e 


ar^yjriiiiriftf Ti.va»nf 



19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 

M, (M 0 ) 

Figure 3. Fraction of runaway stars that escape the supercavity to produce observable SNRs 
This is for R^-i * 5 kpc and includes the gravitational restoring force. 
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Figure 4. Fraction of runaway stars that escape the supercavity to produce observable SNRs. 
This is for * 10 kpc including the gravitational restoring force. 
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TABU 1 


A. The Molecular Cloud 81% (n - 10* cm" 3 , T * 10 2 K) 


t(yrT 


V s (ka a’ 1 ) 


V k) 


Re (PC) 


Remarks 


0 

5000 

3.5x10" 

0 

Explosion occurs 

23 

5000 

3.5x10® 

0.12 

Free expansion ends 

112 

1200 

2xl0 7 

0.3 

Adiabatic phase ends 

9x10* 

v; 5 

2 . 5x1 0 3 

1.9 

Shell stalls 

The Hot Cavity SNR (n - 

10~ 2 cm“ 3 , T 'v. 5 

x 10 5 K) 


t(yr) 

V(km s" 1 ) 

T„(k) 

R-(dc) 

Remarks 

3 a a - 


0 

5000 

3.5x10 

0 

Explosion occurs 

3.7xl0 3 

5000 

3.5x10® 

19 

Free expansion ends 

1.4xl0 5 

300 

1 . 3x10® 

105 

SN shock encounters 





moving bubble shell 





and it quickly gets 





decelerated to the 

.3 

The Classical SNR (n ~ 1 cm 9 

T 10 2 - 

10* K) 

bubble velocity (21 km s 3 ) 

t(yr) 

V(km a -1 ) 

T„(k) 

R„(pc) 

Remarks 

5 9 O 

0 

5000 

3.5x10® 

0 

Explosion occurs 

6.2xl0 2 

5000 

3.5x10® 

3.1 

Free expansion ends 

2.9x10* 

265 

10 6 

19.6 

Adiabatic phase ends 

3.6xl0 6 

5 

2.5xl0 3 

73.6 

Shell stalls 


29 



TABLE 2 


Primary and Secondary Masses of the Binary System 
and the Evolutionary Timescale Differences 


M, (M )* 
1 © 

"2 <M / 

AT 2 i<yr> 

25.85 

19 

2.8xl0 6 

24.49 

18 

2.9xl0 6 

23.13 

17 

3.7xl0 6 

21.77 

16 

4.0xl0 6 

20.4 

15 

4.2xl0 6 

19.05 

14 

5.1xl0 6 

17.69 

13 

5.4xl0 6 

16.33 

12 

5.6x10^ 

15.0 

11 

7.9xl0 6 

13.6 

10 

9.9xl0 6 

12.24 

9 

1.49xl0 7 

10.88 

8 

1.86xl0 7 

9.53 

7 

2.43xl0 7 

8.16 

6 

3.03xl0 ? 

6.8 

5 

4.25xl0 7 

5.44 

4 

7.69xl0 7 


* Average mass of the primary for the tabulated mass of the runaway secondary. 
+ Mass of the runaway secondary. 

** Difference in the total evolutionary times, & T 21 ' T 2~ T 1’ where ^ is 

the evolutionary time of the primary and is the total evolutionary time 

of the runaway secondary. 
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TABLE 3 


A Lower Limit on the Percentage of She 
That Produce Observable SNR*** 


SN Progenitor 
Maes Qf$) 

17-70 

17-16 

16-15 

15-14 

14-13 

13-12 

12-11 

11-10 

10-9 

9-8 


Relative Fraction 
of Total SN 
Progenitors 

17.9 

2.7 

3.3 

4.1 


5.2 

6.7 

8.7 
11.7 

16.3 

23.4 


Percentage of $Ns + 
Producing Observable SNRe 


*C«1 - 5kpc 

Rc.i " 10kpc 

*Cal - 201 

0.0 

0.0 


0.9 

0.0 

% 

1.1 

0.0 

pH 

1.4 

0.0 

f . 

1.6 

0.4* 

ist 

1.8 

1.3 

0 S 

H 

2.0 

2.4 


2.0 

3.0 

6* 

2.1 

3.2 

° ■§ 

2.1 

3.2 



Runaway Total 100. 0Z 

Additional Contribution 
by Slowly Moving Stars 

^Total Allowing for 30% ISM 
as Supercavities 


15. 0Z 

13. 5Z 

O.OZ 

4.3Z 

O.OZ 

O.OZ 



O.OZ 

14. 3Z 

9.0Z 

O.OZ 


*The mass bin is only 13.0 

-a/H 

+A density law n - n Q e 

to 13.5 

is assumed and with 

following values: 

Rg^Ocpc) 

n o (cm‘ 3 ) 

H(pc) 

5 

3 

70 

10 

1 

150 

20 

0.1 

500 


4+These total percentages are computed by assuming that 1/3 of the galactic plane 
is occupied by supercavities, i.e. it is computed by multiplying the Runaway 
Total by 2/3, but all slow moving stars contribute. 

^Assumptions: 8M Q is the lower limit on SN progenitors. 

The gravitational restoring force is negligible. 


V - 50 km s 
P 
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TABLE 4 + 


Mass Bin Percentage of SH 

(Mq) Producing Observable SNR 


R_ . 
Gal 

• 51cpc 

R Cal ’ 10k P C 

19-18 

0 

0 

18-17 

0 

0 

17-16 

0.12 

0 

16-15 

0.23 

0 

15-14 

0.29 

0 

14-13 

0.36 

0.36 

13-12 

0.46 

0.46 

12-11 

0.61 

0.61 

11-10 

0.72 

0.67 

10-9 

1.13 

0.72 

9-8 

1.63 

0.89 

8-7 

2.45 

1.35 

7-6 

3.85 

2.67 

6-5 

6.8 

6.8 

5-4 

12.39 

7.22 

Runaway Total 

31.05% 

21.75 

Slowly Moving Stars 

9.13% 

8.6% 

Total"*"*’ 

29.8% 

23,10% 


"*The same density law as in Table 3 is assumed. 

4 * 4 * 

This total is obtained by multiplying row 16 by 2/3, 

(i.e. we assume as before that 1/3 of all ISM is occupied 
by cavities) and adding it to row 17. 
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